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Abstract: The solution to the Schwinger Dyson equation that describes the summation over Pomeron 
C") loop diagrams is derived. The solution is a closed expression which splits into two parts. The first leads 

directly to the renormalization of the BFKL Pomeron, and the second contribution is equivalent to non 
qs^ interacting Pomerons with renormalized vertices. Thus a closed expression is derived for the sum over 

Pomeron loop diagrams in the perturbative QCD approach, which preserves unitarity. 
>' 



^ Keywords: I BFKL Pomeron, Triple Pomeron vertex, summing Pomeron loops, QCD, Schwinger Dyson 



Cd |equat 



ion 



*Email: jeremy.miller@ist.utl.pt; miller@physics.org 



Contents 

|]. Introduction [I] 

|2|. The bare Pomeron scattering amplitude |^ 

|||. The Schwinger-Dyson equation |5| 

The non interacting Pomeron solution 10 

§. Results [L3 

H. Conclusions and discussion Il3 



1. Introduction 

The goal of this paper is to derive a solution to the 
Schwinger-Dyson equation, which provides the full sum over 
Pomeron loop diagrams. The motivation for pursuing the 
summation over Pomeron loop diagrams, is the large contri- 
bution of Pomeron loops to the diffractive scattering ampli- 
tude in short distance interactions. Hence a reliable calcu- 
lation of the scattering amplitude demands the summation 
over Pomeron loop diagrams to be taken into account. The 
bare scattering amplitude from the t-channel exchange of a 
single Pomeron, and also Pomeron loop diagrams grow with 
energy. Unitarity is only restored by replacing the Pomeron 
Green's function with the sum over the full set of loops. 

The BFKL Pomeron shown in Fig. [I] is the t-channel exchange of a pair of vertical gluons, interacting 
through horizontal gluons which form the "rungs of the ladder" structure. The BFKL Pomeron is the sum 
over all ladder diagrams of this type with n rungs of the ladder. The vertical gluons are themselves a 
superposition of the sum over n rung ladder diagrams, and so on. This leads to the scattering amplitude 
A (s, t) which was derived in refs. [1-5] in the leading log approximation to be proportional to 




gluon BFKL Pomeron 

proton 

Figure 1: The BFKL Pomeron 
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A (s, t) oc V -r (a G (t) In ( S /s )) n = s QG(i) (1.1) 
n! 

n=0 

where a G (i) is the Regge trajectory. In this way the sum over ladder diagrams is achieved by replacing 
the two interacting vertical gluons with a "reggeon" which behaves as s aG W at high energy. According to 
the optical theorem the total cross section behaves as a tot oc s a<3 *- ^ _1 . Experimentally it is known that the 
total cross section rises slowly with s which means olq (0) > 1. Pomeranchuk first commented [6, 7] that 
this behavior is matched by the theoretical prediction that ac (0) > 1 when the t-channel exchange carries 
zero quantum numbers, including zero charge and color flow. Such particles with quantum numbers of the 
vacuum exist in QCD for bound gluon states. This kind of bound state trajectory is called the Pomeron 
named after Pomeranchuk, which is the double t-channel gluon exchange shown in Fig. ffl. The evolution 
of the vertical t-channel gluons to the sum over ladder diagrams is called the BFKL Pomeron which is 
described by the BFKL equation [8, 9]. 

For hard collisions, Pomeron loop corrections contribute substantially to the scattering amplitude, and 
a reliable estimate of the scattering amplitude requires the sum over Pomeron loop diagrams. This is a 
difficult problem, and for a long time the generally accepted method for estimating the sum over Pomeron 
loop diagrams was the Mueller, Patel, Salam and Iancu (MPSI) approach, (see refs. [10-16]). A. Mueller 
[15] and Levin, et. al. [16] first commented that at high energy Pomeron loop diagrams should reduce to 
independent Pomeron exchanges with complicated non factorized vertices. In refs. [17, 18] it was shown 
that this problem can be solved theoretically in perturbative QCD for a specific type of diagrams. In this 
approach the special class of symmetric diagrams of the type shown in Fig. ^ were calculated using an 
iterative technique, based on the observation that Fig. ||] (b) is generated from Fig. || (a) when each branch 
of the loop gives birth to a secondary loop leading to two "second generation" loops. Likewise Fig. || (c) 
arises when when each of the two "second generation" loops in Fig. || (b) gives birth to two loops which 
leads to four "third generation" of loops. Hence the diagrams in Fig. ^ are called the N = 1, N = 2 
and N = 3 generation loop diagrams, and continuing in this way one can generate the full spectrum of 
symmetric Pomeron loop diagrams, with N generations of loops. 

The formula derived in ref. [17, 18] for the sum over this class of symmetric diagrams was of the 
type Y2n^ (-^0> w h ere A (N) is the diagram with N generations of loops. This formula suffered from the 
difficulty that it diverges with energy. In this paper, this problem of unitarity violation is resolved by 
taking instead the sum ^ n A (n) where A (n) is the symmetric diagram with n Pomeron loops. In this 
way a closed analytic expression is obtained, which is equivalent to the sum over diagrams with In non 
interacting Pomerons. 

However this does not complete the sum over Pomeron loops, since the class of diagrams with succes- 
sive loops shown in Fig. ^ should also be included. It was first suggested by M. Braun in refs. [19, 20] that 
the Schwinger-Dyson equation automatically generates the sum over the full spectrum of Pomeron loop 
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diagrams, which includes both types of diagrams shown in Fig. || and Fig. ||[ 



This paper is organized in the following way. Firstly in section [2] the scattering amplitude arising from 
a t-channel bare Pomeron shown in Fig. [l] is derived, for the sake of completeness. Next a solution to the 
Schwinger-Dyson equation is presented which splits into parts. The first part discussed in section |3| leads 
directly to a simple expression for the renormalized Pomeron intercept, which arises from summing over 
the class of diagrams in Fig. The second part of the solution derived in section || is equivalent to the sum 
over non interacting Pomeron diagrams, which is derived from the sum over the type of diagrams in Fig. ^. 
Intuitively this can be seen from an observation of Fig. ||, where at high energy taking the branches of the 
loop in Fig. |3| (a) outside leads to 2 non interacting Pomerons. Likewise for the 2 small loops in Fig. || (b) 
that have not given birth to any more smaller loops, taking the branches of the loops outside leads to 4 
non interacting Pomerons. In section || the main results of the paper are presented, and in section ^ the 
conclusions of this paper are discussed. 



2. The bare Pomeron scattering amplitude 

In this section the diffractive scattering amplitude arising from a t-channel bare Pomeron shown in Fig. |l], 
is derived. Although the expression is well known and has been calculated in refs. [16-18, 21-25], it has 
been included for the sake of completeness. The solution to the BFKL equation provides the trajectory 
u(n,u) (neZ;i/£R) for the BFKL Pomeron as [8, 9]; 



u(n,v) = a 8 U(l)-^i(>n^ + iujy, a s = ^ (2.1) 

where ip (x) = dlnT (x) jdx is the di-gamma function, n represents the energy levels of the BFKL 
Pomeron, and v is a continuous variable which one integrates over when calculating Feynman diagrams. 
The BFKL eigenfunction falls sharply with increasing n and is only positive at high energy when n = 0. 
Hence throughout this paper which is focussed on high energy scattering, n = is assumed and the 
argument n is suppressed. Hence the BFKL Pomeron trajectory which is the sum over ladder diagrams 
of the type shown in Fig. |l| is described by the regge behavior = e w ( u ' v . The scattering amplitude of 
Fig. [I] is given by the expression; 



A (0) (y s % r |Fig.i = -^/ dvh(y)g{v)e^E v E'_ v (2.2) 

h{v) = ±y. g(u) = - E v =(^-\ [-P^-] (2.3) 

K> ^ {J 16(l/4 + ^) 2 ' \r w r 20 J \r* 10 r*J { ' 

h{v) is the integration measure which preserves conformal invariance [19, 20], g (u) is the Pomeron 
propagator in the conformal basis [19, 20] and E„ is the coupling of the BFKL Pomeron to the QCD color 
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dipole [19-21], in the dipole approach to proton proton scattering. Here r%2 = r\ — r<i is the transverse size 
of the dipole and t\q = t\ — tq where ro is the center of mass coordinate of the dipole. The observation that 
the BFKL eigenfunction Eq. ( |2.lD has a saddle point v = means that one can expand the exponential in 
Eq. flO) as 



(v) = to (0) - -v 2 u" (0) + O {u 2 ) ; w(0) = 4a s ln2; uj" (0) = 28<5 S C (3) 



(2.4) 



where C, (3) = 1.202 is the Riemann zeta function. Using this expansion the integration in Eq. ( |2.2| ) is 
evaluated by the steepest descent method which yields the result [17, 18]; 



4(0) {y,Sy H \Fig. |) 



a 2 (2vr) 



1/2 



2tt 2 N? 



(u;"(0)yf 2 



(2.5) 



3. The Schwinger-Dyson equation 





i+i 
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Figure 2: Diagram (a) shows n Pomeron loops in se- 
ries. Diagram (b) shows a series of n Pomeron self mass 
interactions, where the bold lines represent a superposi- 
tion of n loops in series. 



Eq. (|2.5[) is the scattering amplitude arising from 
the exchange of a bare Pomeron. According to the 
Schwinger-Dyson equation [19] , the bare Pomeron prop- 
agator should be replaced by the Green's function 
which is found by summing over the class of diagrams 
of Fig. ||. This Green's function is found by first sum- 
ming over all diagrams with n consecutive loops in se- 
ries shown in Fig. Q (a), from n = to infinity. Next, 
the Pomeron lines themselves in Fig. |2| (a) should 
be replaced with the sum over loops in series which 
leads to Fig. § (b). The self mass terms in Fig. § (b) 
are "loops made of a superposition of loops" . Con- 
tinuing in the same way to introduce more gener- 
ations of loop series, one derives the Green's func- 
tion which is the sum over all Pomeron loop dia- 
grams. The process which has been explained here 
in words, is described by the Schwinger Dyson equa- 
tion introduced by M. Braun in ref. [19, 20]. This 
expresses the full Pomeron Green function G (v,y) in 
v, y-representation, as the sum over the class of dia- 
grams in Fig. |2| as; 
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'•y ryi 

G{v,y) = g(v) - g(u) I dy 1 / dy 2 m{y,y l2 )G [y,y 2 ) (3.1) 



o J o 



where m fa y\ 2 ) is the Pomeron self mass {y\ 2 =y\ — y 2 ), which sits between the rapidity values y\ for 
the upper rapidity limit, and y 2 for the lower rapidity limit (see Fig. |2| (a)), and is given by the expression; 



m 



-i poo poo 

fal/12) =-/ du 1 h(u 1 )G(u 1 ,y 12 ) du 2 h (y 2 ) G fa, y 12 ) \F fa u u v 2 ) | 2 e (^M+^)-^))yi2 (3 . 2) 

^ J-oo J — oo 



The pre-factor of 1/16 in Eq. ( |3.2| ) divides by the order of the symmetry group, so that identical dia- 
grams are only counted once (see ref. [17] for a full explanation). T fa v\, u 2 ) is the triple Pomeron vertex 
for the splitting of the Pomeron with the conformal variable v into two daughter Pomerons with conformal 
variables v\ and v 2 which form the branches of the loop. The splitting vertex is the complex conjugate of 
the merging vertex, so the squared absolute value of the vertex in Eq. ( |3.2| ) is the product of the splitting 
and the re-merging vertex, forming the loop. 



The Schwinger Dyson Eq. ( |3.1| ) forms an iterative sum, which after expanding takes the following form; 

00 n fy l 2 +1 ry\ 

G fa y) = g fa £ (-1)™ g n (u) ]J dy\ dy\ m (u, y\ 2 ) ; y 2 n+1 = y; (3.3a) 

n=0 i=l J ° J ° 

1 f'OO f'QO 

m{u,y\ 2 ) =- d^hWGfayiJ dv 2 h (y 2 ) G (u 2 , y\ 2 ) \T {v, u u v 2 ) | 2 e («W^H("))iS a ( 3 .3b) 

J —oo J —oo 

where m (v, y\ 2 ) is the Pomeron self mass which sits between the rapidity values y\ and y l 2 in Fig. || (b), 
and y\ 2 = y\ — y 2 - The G (^1,2/12) and G (v2,y\ 2 ) in Eq. ( 3.3b| ) are the same Pomeron propagators given 



by the infinite Schwinger-Dyson sum of Eq. ( |3.3a| ). From this it becomes clear how plugging Eq. ( p. 3b ) into 



Eq. ( 3.3a| ) generates the non-closed sum over the never ending spectrum of Pomeron loop diagrams. The 



integration limits are due to the upper rapidity value y\ in Fig. g which cannot exceed the lower rapidity 
value y % 2 +l of the next self mass interaction above it. The strategy for deriving a closed expression which 
is a solution to Eq. ( |3.3a ) is the following. 



First consider the sum over diagrams of the type shown in Fig. || (a), which contains n Pomeron loops 
in series. The sum over all such diagrams from n = to infinity is described by the Schwinger Dyson 
Eq. (|3.3a ) by replacing m (v, y l l2 ) with the Pomeron loop tuq (v, y l l2 ) which yields; 



G (^y|Fig.g(a)) =g(v)J2(-l) n g n (u)H / dy\ dy\m {v,y\ 2 ) (3.4a) 

ZTn JO JO 



n=0 i=l 

00 



1 poo poo 

mo (v, y\ 2 ) = -t / dvth {yx) g fa) / du 2 h fa) g fa) |r (y v u v 2 ) | 2 e M-i)+<^)-"WM 2 (3.4b) 
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where mo (^,2/12) is found from Eq. ( |3,3b[ ) by replacing the renormalized propagators G (z/i, yf 2 ) anc ^ 
^(^2,2/12)) w hh the bare propagators g(vi) and 3(^2)- From the Korchemsky expression for the triple 
Pomeron vertex [26], it was found in ref. [17, 18] that there are two main contributions to the triple 
Pomeron vertex which are given by the following asymptotes; 



, f 1 \ l2 (47r) 6 at 1 

ir (i» = = o,„* = oj l 2 = Sv^T74T^ (3 ' 5a) 

X (i/) = ReU(l)-^Q + ii/)] (3.5c) 



Eq. ( ]3.5a| ) leads to the contribution to the Pomeron loop amplitude of Eq. ( |3.4b| ) which is equivalent to 2 
non interacting Pomerons, with renormalized Pomeron vertices. This particular contribution to mo (u, 2/12) 
gives a vanishing result for n > 1 in Eq. ( 3.4a ). This makes sense since Eq. ( 3.4a ) describes the sum over 
n Pomeron loops in series shown in Fig. ^ (a). Therefore the phenomena where the branches of the loop 
span the entire rapidity gap between the projectile and target to become 2 independent partons, cannot 
occur with more than 1 loop in series. With this in mind the contribution to the ver tex of Eq. flpD 

is 

postponed until the non interacting Pomeron solution is discussed later on in section ||. For now, inserting 
the asymptote of Eq. ( 3.5b| ) into Eq. ( |3.4b ), and then introducing the definitions given in Eq. ( p. 3] ) (where 
g (u) can be cast as 1/16 (1/2 + iv) 2 (1/2 — iv) 2 ); 




The v\ , V2 integrals are solved by closing the contour over the upper half plane and summing over the 

residues at iv\, ii>2 = 1/2, taking into account the poles which stem from u (^1,2) — — — -a s / (1/2 — iv\ 2) 
such that Eq. (|3.6|) becomes 
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mo \y 



yj 2 ) 



7T 

y 



i 



c 



dv2 v\ (1 — %V\ — Wq) 



xe 



C (i + ii/i) 2 (i-ii/i) 2 \a S dy\ 2 , 



di^i/ 2 (1 - m - iv 2 ) c (-a s{ H2~i V2 r x -^M)vU (3 7 ) 

C {\ + iv 2 f {\-iv2f 

It is instructive to change the v\ integration variable to u = a s / (1/2 — iv\) + a s / (1/2 — iv 2 \ such that 
the jacobian cancels the remaining singularity which stems from (1/2 — iui)~ 2 in Eq. Q3.?1 ). Integrating over 
u yields the derivative of the Dirac delta function 2iri5^ (y\ — y 2 )- After taking the residue at %v\ = 1/2 
the v 2 integral reduces to § c dv 2 v 2 / (5 + ^2) 2 Q — w%) 3 = 27ri/4. Over all, after re-arranging the order of 
derivatives so that there are no derivatives of the delta function in the integrand, one finds for the Pomeron 
loop amplitude the following expression; 



m (u, y\ 2 ) = a Q + A ' x (u) u? (y) e'^M (y\ - y\) ■ 



a = a '\ 1 " Jp ) ( 3 - 



Finally inserting Eq. ( p.8[ ) into Eq. (3.4a), the integrations over the rapidity variables are trivially 
solved thanks to the Dirac delta function. Thus one derives the following expression for the sum over the 
class of diagrams of Fig. [2] (a) with n consecutive loops ; 



00 (—'\\ n 

Go (u, y|Fig. I (a)) = g (y) J] { -^- (g (u) a (l/4 + v 2 f x („) u? (y) y 



n=0 



• ry(//)l ' xp ( ~TE (l + i/2 ) x{y)uj2 (,/)// ' 



Substituting for the bare propagator g (v) that appears in Eq. ( |2.2| ), the renormalized one derived in 
Eq. (|3.9p, leads to the scattering amplitude which is equivalent to the replacement; 



w („) - 2 („) = u (u)+t {v) ; C 0") = (~ + X (1/) ^ 2 M (3.10) 

Therefore the sum over the class of diagrams in Fig. |2| (a) leads directly to the renormalized Pomeron 
intercept given by Eq. ( gig ). Phrased differently, the sum over the class of diagrams of Fig. y (a) is found 
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by replacing uj (is) in Eq. ( |2.2| ) with Eq. fl3.10| ). The bold lines in Fig. [2] (b), represent the sum over the 
class of diagrams in Fig. |(a). As such Fig. |(b) is a string of loops, where the bold lines which form the 
loops are themselves a string of loops. The bold lines are reggeons, with the renormalized intercept uj (z/) 
derived in Eq. ( 3.10| ). With this in mind, the sum over the class of diagrams of Fig. |2| (b) is derived from 
Eq. ( pOJ ), by replacing the intercepts uj (vi) and uj {v 2 ) with uj and uj (v 2 ) found in Eq. ( 3.1 0| ) , which 
yields; 



00 n rvi +1 rv\ 

G (u, y|Fig. 1(b)) = 9 (u) £ (-1)" g n (*/) H dy\ dy\ m (u, y\ 2 ) ; (3.11a 

n=0 7=1 J ° J ° 

|2 „(S(i^)+w(w2)-w(i/))j/: 



j /■oo 

16 



V12) = T7 / ^1/1 M 5 (fl) / ^ M <7 M |r (!/, zvi, */ 2 ) \ z eS"KWW"WM2 (3.11b) 



Following the same arguments used above, the integrations in Eq. (3.11b) are solved by closing the 
v%, v 2 contours over the upper half plane and summing over the residues at iv\ y iv 2 = 1/2. After introducing 
the definitions given in Eq. ( |2.3| ) and taking into account the singularities that stem from the triple Pomeron 
vertex of Eq. ( 3.5b|) , and the asymptotes; 



uj{y) 



U/-+1/2 — a s 



CM 



16 (1 



(3.12) 



Eq. (3.11b) becomes; 



m \ v 



,2/12) 



7rV _J_ 

TV 1 " m 



2 1 exp , 1 

c {\+ivx) {\-ivi) \ 2 " 



a s ,• a a 
2/12 + 



C (|+^2) 2 (|-^ 2 ) 



3 (w(iAa)-a)(i/))i/J 2 



16 (|--l) 2 



2/12 



7T 



-1 d 



dv 2 v 2 (1 — iv\ — ivi 



G (|+^2) 2 (|-^2) 



3 (w(i^)-a)(i/))i/J a 



a i d 



exp — y 12 



-a. 



16 a ' lz d 2 y\ 2 J ^ 6XP I i — ii/i 



2/12 



(3.13) 
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Changing the integration variable tow = a s / (1/2 — ivi)+oj (i^) and integrating over w yields the Dirac 
delta function 2iri5 — y^) , acted on by the derivatives with respect to y\ 2 which appear in the integrand. 
After taking the residue at iv\ = 1/2 the U2 integral reduces to § c dv2 v \l (5 + ^2) 2 {\ ~ ^"if = 2vrz/4. 
After rearranging the order of derivatives so that there is no derivative of the Dirac Delta function in the 
final expression; 



m 



M2) = a ( L A +v 2 ) xWw 2 M e-^{exp 



2 a/1 - 4ip (u)ip 2 (y) Id 



<f>(u)}6(y{-yi) (3.14a) 
(3.14b) 



Finally inserting Eq. ( 3.14a| ) into Eq. (3.11a) and evaluating all the integrations over the rapidity 
variables leads to the following result for the sum over the class of diagrams shown in Fig. (b); 



afl 2 



G (u, y|Fig. 1 (b)) = g (u) exp { -- I - + 1/ ) X iy) uj z (u) (1 + (u)) y } (3.15) 



Using exactly the same arguments which led to Eq. ( |3.10| ), the sum over the class of diagrams with 
consecutive self mass terms shown in Fig. |2| (b), is achieved by replacing u> (y) that appears in Eq. (|2.2| ) 
with the following renormalized Pomeron intercept; 



n 



u(y)-^u (y) =uj(v) + Q (v) (1 + (v)) 



(3.16) 



If the same treatment is repeated and another deeper level of loops are introduced, the Pomeron 
intercepts ui (ui) and oj (^2) in Eq. (|3.4b| ) are replaced by Eq. ( |3.16|) . The u±, V2 integrals are solved in the 
same way by closing the contour over the upper half plane and summing over the residues at iv\,iv2 = 1/2. 
Fortunately, ^(^1,2) vanishes as iv\^ — > 1/2, so it gives no contribution to the sum over residues. This 
means that the steps from Eq. ( 3.11| ) to Eq. ( |3. 15 ) are identical, and the same expression of Eq. ( 3.16] ) for 
the renormalized Pomeron intercept is derived. Therefore every time a deeper level of loops is introduced, 
one always arrives at Eq. ( 3.16| ) for the renormalization of the Pomeron intercept. Therefore the following 
expression derived in Eq. ( 3.15[ ) 



(3.17) 



is the solution to the Schwinger Dyson equation Eq. (3.3a) which describes the sum over the class of 
Pomeron loop diagrams in Fig. [| The solution is equivalent to the replacement of the Pomeron intercept 
oj (f ), with the renormalized Pomeron intercept oj (y) derived in Eq. ( |3.16 ). 
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4. The non interacting Pomeron solution 




Figure 3: The special class of symmetric Pomeron Figure 4: The N generation diagram which stems 

loop diagrams taken into account in the sum over from the simple loop giving birth to two sets of N — 

Pomeron loops, (a) is the diagram with N = 1 gen- 1 generations of loops, is equivalent to the diagram 

eration of loops, (b) has N = 2 generations of loops of 2 N non interacting Pomerons with renormalized 

and (c) has N = 3 generations of loops. Pomeron vertices. 



The next part of the discussions is focussed on the non interacting Pomeron solution, which stems from 
diagrams of the type shown in Fig. [3|. Loops which become non interacting Pomerons do not contribute to 
consecutive loops in Fig. [|| (a), because the phenomena where the branches of the loop stretch and span 
the whole rapidity gap between the projectile and target, can't happen for more than one loop in series. 
Therefore this solution cannot be treated in the context of the Schwinger Dyson Eq. ( j3.3a|) for n > 1, and 
requires a separate approach. Using the same conventions, the scattering amplitude of Fig. ||| (a) is given 
by the expression; 



IT 2 f°° 

A (1) (y|Fig. 1 (a)) = -f / duh (u) g 2 (u) e^m (1) (u, y) E V E'_ V (4.1a) 

f'OO poo 

m {1) (v,y) = —l dvih (vi) g (vi) I dv 2 h (u 2 ) g (u 2 ) \T u 2 ) | 2 (4.1b) 

ro / — / — r>r . 



y pyi 

(u)(v\)+u)(v2)-u(v))y\2 



dyi / dy 2 e ( 



o 



where y± 2 = y±—y 2 is the rapidity gap which the loop fills (see Fig. || (a)). Now inserting the asymptote 
of Eq. ( 3.5a| ) for the triple Pomeron vertex for the region where is\, v 2 are close to zero, one can substitute 
for the BFKL eigenfunctions uj (yi) and uj [y 2 ) the expansion of Eq. ( |2.4D and integrate over v\ and v 2 using 
the method of steepest descents, which gives the expression [17, 18]; 



m 



(i) 



y) 



(1/4 + 



3 (2a)(0)-w(i/))yi2 



dyi / dy 2 - 



2 10 a! 



Vl2 



(4.2) 



Inserting Eq. (4^) back into Eq. ( |4.1a| ), the v integration can be solved by closing the contour over 
the upper half plane and summing over the residues at iv = 1/2, taking into account the singularity which 
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stems from uj (y) given in Eq. ( |3.12| ). This leads to the following contribution to the scattering amplitude 
[17, 18]; 



Eq. (4.3) is equivalent to 2 non interacting Pomerons, with renormalized Pomeron vertices. As ex- 
plained above, Fig. ||| (b) stems from Fig. [| (a) when each branch of the loop gives birth to a secondary loop 
leading to the two "second generation" of loops in Fig. || (b). In the same way when the second generation 
loops in Fig. || (b) each give birth to two loops, this leads to 4 "third generation" of loops in Fig. || (c). 
Continuing with this evolution, the entire spectrum of symmetric N generation diagrams can be generated 
for all N . The scattering amplitude with N generations of loops shown in Fig. ^| is the generalization of 
Eq. ( gig ), namely [17, 18]; 



rv 2 f 00 

A {N) (y|Fig. 1) = -f duh{u)g\u)e^ym {N) {v,y)E v E'_ v (4.4) 

^ J — oo 

where mny) { v i v) ls the contribution of the N generations of loops in Fig. ||. In refs. [17, 18] a detailed 
explanation of how to calculate rrifm (u, y) was given. This is based on the observation that Fig. || is 
equivalent to the simple loop diagram of Fig. || (a) when each branch in the loop gives birth to a set of 
iV — 1 generations of loops. This means that to write the expression for mny) (y, y), all that is needed is to 
modify the propagators for the branches of the loop in the expression of Eq. ( |4.1b[ ) as; 



#0i) -»■ 9 (^i) "i(jv-i) {"uvvd 9{vx) ; gfo) -»■ g iy 2 ) m^ N _^ (u 2 , 2/12) g (^2) (4.5) 



After implementing Eq. ( |4.5| ) ha Eq. ( 4.1b| ), one arrives at the following amplitude for the set of N 
generations of loops; 



16 



00 



m, N) (v\y,5y H ) = -rx I du x h (fi) g (fi) I du 2 h (u 2 ) g (u 2 ) \T {u\v\, v 2 ) \ 2 (4.6) 
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Eq. (4.6) forms an iterative expression. Using the technique of proof by induction, the following formula 
derived in refs. [17, 18] for Eq. (|4.6j ) can be proved; 



2[JV " 1] n 



b> 

2 10 af 



yi 



iV>[w"(0)] 3 ' 



t + v j x (» J dy 1 J dy 2 y 12 



2 11 V JV? 



2 [iv-i]_ 1 



1 d \ 3 e M%l2 



2 [JV-1] 



y? 2 



(4.7a) 



(4.7b) 



Finally after inserting Eq. Q4.7D into Eq. (4.4), one finds the following scattering amplitude for the 
diagram of Fig. with N generations of loops; 



A 



(TV) 



|Fig. § 



a s {bb ) 2 [n-i] 
T, — V 



7T 



-1 d\ 3 e 2 ^y 

a s dyj y 3 



2 [N-1] 



(4.8 



Eq. (4^8) is equivalent to 2^ non interacting Pomerons, with renormalized Pomeron vertices shown 
pictorially in Fig. 0. This follows from the observation that Eq. fl4.8|) can be recast in the form; 



A(N) (y|Fig. 



(4.9) 



where the coefficient k (n) contains the set of renormalized Pomeron vertices. Eq. (4.8) describes the 
scattering amplitude with 2 N ~ 1 loops, which is equivalent to 2^ non interacting Pomerons. This formula 
can be generalized to the scattering amplitude of the symmetric diagram which contains n loops, or 
equivalently 2n non interacting Pomerons, namely; 



A (n) (y) 



(bb'y 



2 9 iV c 2 vr 



-1 d\ 3 e 2uj(0)y \ 
a s dyj y 3 i 



(4.10) 



The sum over the complete set of symmetric Pomeron loop diagrams, is achieved by evaluating the 
sum Y^=o A {n) (y I Eq. (|4.10D ). In the outcome formula, the intercept ui (0) should be replaced with S3 (0) 
of Eq. ( |3. 16 ) . This leads to the sum over diagrams with independent Pomeron exchanges, where the 
Pomerons are replaced by the superposition of loop series shown in Fig. ||, described by the Schwinger 
Dyson equation. In this formalism the full sum over symmetric Pomeron loop diagrams is given by; 



n=l 



2 9 iV c 2 ^ 



. ( 1 d \ 3 e 2 ^y 

a s dyj y 3 



(4.11) 
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^ a s dyj y 3 



- 12 - 



5. Results 



1.0 h 



LA (ylEq.4.11 ) 
(n) 



Rapidity y 



The energy dependence of the 
symmetric Pomeron loop diagrams 



Figure 5 

sum over 

EZiAn) (y|Eq. fluij)) 



The energy dependence of the sum over Pomeron loop 
diagrams is shown in Fig. [|. The curve approaches the 
black disk limit, so that unitarity is preserved. If the bare 



Pomeron intercept uj (0) is used in Eq. (4.11) instead of the 
renormalized one, the graph of Fig. [5] is unaffected. This 
indicates that the dominant contribution to the sum over 
Pomeron loop diagrams comes from the diagrams which are 
equivalent to non interacting Pomerons, with renormalized 
Pomeron vertices. The class of diagrams of Fig. ^ leading to 
the renormalized Pomeron intercept give a negligible con- 
tribution in comparison. 



The formula of Eq. ( |4 . 1 1| ) requires explanation. It is tempting to think that Eq. ( 4.11[ ) is the scattering 
amplitude, however a closer look reveals that this formula is still just the sum over a special class of loop 
diagrams, as will now be explained. 



6. Conclusions and discussion 

The following discussion is a summary of the formalism adopted in this paper, for the summation of 
Pomeron loops. There are two distinct ways of summing over Pomeron loops, and both methods must be 
taken into account, namely; 



1. The sum over the class of loops in series of Fig. ||, using the Schwinger Dyson Eq. (3.3). The asymptote 



for the vertex of Eq. ( 3.5b| ) contributes to the Schwinger Dyson sum, whereas the asymptote of 



Eq. ( 3.5a| ) gives a vanishing contribution for n > 1. The Pomeron loop summation generated by the 



Schwinger Dyson sum, provides the renormalized Pomeron intercept derived in Eq. ( |3.16| ). 

The sum over the symmetric class of loops in Fig. ^, using the vertex of Eq. ( |3.5a| ). This leads 
to the sum over even numbers of non interacting Pomerons, with renormalized Pomeron vertices. 
The symmetric nature of the loops in Fig. [3|, leads to even numbers of independent Pomerons. For 
example taking all loop branches outside in Fig. ||| (b), leaves 4 independent Pomerons. 



The two asymptotic expressions for the triple Pomeron vertex in Eq. ( |3.5| ), lead to the above two 
entirely different types of loop summation. It should be stressed, that the above two treatments do not 
lead to the same result, and the scattering amplitude requires taking into account both methods. The non 
interacting Pomeron solution (2) , originates from the loops in Fig. ^ stretching in rapidity space until they 
fill up the gap between the projectile and target, and therefore become independent Pomeron exchanges. 
This phenomena does not occur for the loops in series shown in Fig. B, since the loop cannot become non 
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interacting Pomerons, when there is more than one loop in series. The only class of diagram which can 
yield non interacting Pomerons, is the special class of loops shown in Fig. || (a), and the only asymptote 
for the vertex which can yield non interacting Pomerons is Eq. ( 3.5a| ). Hence, the non interacting Pomeron 
solution requires a separate treatment from the Schwinger Dyson equation. 




+ 



(b) 




+ 



+ 



(c) 

Replace co(v) 
with the renormalized 
Pomeron intercept 
of Eq. (3.16); 

co(v) — - co(v) 






Bold Pomerons denote 
the Schwinger Dyson 
sum over diagrams in 
in Fig. 2 . 



Figure 6: The process which leads to the formula of Eq. (4.11). Summing over the symmetric class of loops in (a), 
leads to the sum over even numbers of non interacting Pomerons shown in (b) . Then replace the bare Pomerons with 
bold Pomerons, where bold Pomerons label the Schwinger Dyson sum over the class of loops in series shown in Fig. |[ 
The Schwinger Dyson sum renormalizes the Pomeron, by replacing the Pomeron intercept with the renormalized 



Pomeron intercept lo (v) of Eq. ( 3. 16 ). This leads to the sum over even numbers of independent bold Pomerons in 
(c), that are renormalized in the framework of the Schwinger Dyson equation. 



The result of Eq. ( 4.11[ ) includes both types of summation over Pomeron loops. Fig. |6| shows a picture 
for the derivation of Eq. ( [4.11 ). It is instructive to first consider the class of loops in Fig. ^ (a). As 
explained in section the summation over the class of loops in Fig. || (a) is equivalent to the sum over 
diagrams with an even number of non interacting Pomerons, shown in Fig. || (b). Next, replace the bare 
Pomerons with Pomerons which are renormalized by the Schwinger Dyson sum. This means that for each 
of the non interacting Pomerons in Fig. |6| (b), replace it with the sum over loops in series shown in Fig. 
generated by the Schwinger Dyson equation. This is achieved by replacing the bare Pomeron intercept 
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oj (V), with the renormalized intercept oj (y) derived in Eq. ( 3. 16| ) ■ This leads to the diagram of Fig. || (c), 
where the bold Pomerons label the above described renormalized Pomerons with the intercept oj {y). 



The end result of Fig. ^ (c) is a true description of Eq. ( t4.1l|) , namely the sum over diagrams with an 
even number of independent Pomeron exchanges, which are renormalized in the context of the Schwinger 
Dyson equation. The formula of Eq. ( 4.11 ) does not include the bare scattering amplitude of Fig. |l| given 
by Eq. (|2.5|). The complete scattering amplitude which includes the loop corrections shown in Fig. |6| is 
found by adding to Eq. ( 4.11 ), the bare scattering amplitude of Eq. ([2lj|), This would lead to a divergent 
result which violates unitarity. The only remedy for this problem, is to repeat the procedure shown in 
Fig. ||, for the non-symmetric class of diagrams shown in Fig. ^. Although the diagrams of Fig. [7| have been 
included in the Schwinger Dyson sum, this was performed using the vertex of Eq. ( 3.5b| ). The approach 
here, is instead to use the vertex of Eq. ( |3.5a| ) to sum over the non symmetric diagrams in Fig. [?]. This will 
lead to the sum over odd numbers of non interacting Pomerons, that do not contribute to the Schwinger 
Dyson sum. 



This stems from the non symmetric nature of the loops 
in Fig. [7], where for example taking all branches of the loops 
in Fig. (a) outside, reduces the diagram to 3 independent 
Pomeron exchanges. Following the same strategy described 
above, the independent Pomerons are renormalized by the 
Schwinger Dyson sum, by replacing the Pomeron intercept 
Figure 7: Examples of non symmetric Pomeron with the renormalized intercept found in Eq. (|1§. Over- 
loop diagrams. au this yields the sum over odd numbers of renormalized 

non interacting Pomerons. 




Finally, adding the sum over odd numbers, to the sum over even numbers of renormalized Pomerons 
already derived in Eq. (4.11), leads to the expression which takes the following form; 



A( n ) (y|symmetric + non symmetric diagrams) 



c e 



£(0)2/ 



n=0 



1 + de w{ V)y 



(6.1) 



where c and d contain all the other terms which are part of the scattering amplitude. The key property 
of Eq. (|6.l| ), is that it includes the basic amplitude of Fig. |], and it preserves unitarity generating a similar 
curve to the one in Fig. H Eq. fl6.1|) is the pp elastic scattering amplitude, including the full set of Pomeron 
loop corrections. Unfortunately, we have not yet been able to calculate the class of non symmetric diagrams 
shown in Fig. [7], however this work is in progress. In light of this discussion, the prospects for arriving 
at an expression which preserves unitarity, and includes symmetric and non symmetric loop diagrams, are 
hopeful. 
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(a) 

Figure 8: Diagram (a) shows the lowest order 
triple Pomeron vertex, and diagram (b) shows the 
first order correction to the vertex. 



In the calculations performed in this paper, the di- 
agrams which contribute to the vertex in the framework 
of the Schwinger-Dyson equation, were not taken into ac- 
count. Only the lowest order vertex shown in Fig. |8| (a), 
was included in the above performed calculations. To illus- 
trate one example, the first order correction to the vertex is 
shown in Fig. || (b). The full triple Pomeron vertex which 
includes the complete set of corrections is described by the 
Schwinger-Dyson equation for the vertex; 



T (y,v u v 2 \y x ) = r (0) (v,v x ,v 2 ) 



(6.2a) 



I 2/23 
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where 



oo 



Vv 



■Xt 



dvh (y) g (y) 



(6.2b) 



where T^) is the lowest order vertex shown in Fig. |8] (a), and f is the full vertex which includes 
the complete set of vertex corrections described by the Schwinger-Dyson equation. The Schwinger-Dyson 
Eq. ( |6.2[ ) for the vertex is much more complicated in comparison to the Schwinger-Dyson Eq. ([O]) for 
the Pomeron Green function. The author acknowledges, that the set of corrections to the vertex, are also 
required for the formula for the scattering amplitude which includes all possible corrections. This problem 
is very challenging owing to the complexity of the non-closed equation for the vertex of Eq. ( |6.2D , but 
nevertheless attempts to solve this problem are in progress. The two expressions used for the lowest order 
triple Pomeron vertex in Eq. ( |3.5| ), indicate that the vertex is less than unity. Thus, although the triple 
Pomeron vertex was not taken into account in the framework of the Schwinger-Dyson Eq. ([T^) , since the 
vertex is less than 1, corrections to the vertex are expected to give a small contribution. 

In summary, the main achievements of this article include the following; 

1. A closed solution to the Schwinger-Dyson equation in perturbative QCD, which generates the sum- 
mation over the full set of Pomeron loops, leading to the renormalized Pomeron intercept. 

2. A closed expression for the summation over a special class of Pomeron loop diagrams, equivalent to 
non interacting Pomerons which preserves unitarity. 

Both of these achievements are original, and provide a strong foundation for calculating the scattering 
amplitude in perturbative QCD. The remaining corrections which are required for the scattering amplitude, 
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include the non symmetric loop diagrams of Fig. [7| and the corrections to the vertex in the formalism of the 
Schwinger-Dyson equation. The calculation of these additional required corrections, is the next challenging 
problem to be solved. 
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